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Abstract 


A  method  of  analysis  for  determin¬ 
ing  crack  propagation  loads  on  wood 
beams  with  end  splits  is  presented. 
The  method  is  based  on  (i)  linear 
elastic  orthotropic  fracture 
mechanics  concepts,  (2)  theory  of 
Complex  variables,  and  (3)  least 
squares  boundary  value  collocation 
(BVC).  Using  this  method,  the  critical 
stress  intensity  factor,-Kn£-rTt — 
-Douglas  fii  bedirr5~wtttY-end-spUta  is 
determined, to  be  g,035  pal  -  — 

Also,  a  simple  failure  equation  for 
end  split  beams  is  proposed. 


Acknowledgment 

The  author  Is  indebted  to  C.  B. 
Norris  and  E.  C.  O.  Erickson,  who 
conducted  the  experimental  studies 
reported  herein  In  1951.  They 
reported  their  results  in  an  unpub¬ 
lished  report.  Their  results  have  been 
reevaluated  in  the  fight  of  fracture 
mechanics  theory  In  the  preparation 
of  this  paper. 


United  States 
Department  of 
Agriculture 

Forest  Service 

Forest 

Products 

Laboratory' 

Research 
Paper 
FPL  347 


3fStrength  of  Wood 
1  Beams  y/ith  End  Splits, 


3 


JOSEPH  F/lgURPHV/  Engineer 


,T,c.t  Service..  ■cc'.'".vo'Vi  fapcv^ 


oSilZl/  L 


Jiccui  .jion  for 

1'^  o  miMti 
D2Z  IA3 
U:i  nn.'-Jnced 
Ju  ■  i:‘ ic.it Ion 


■ '  y  £o<t«»s 
• :  L  ■i:«d/or 
-  v-'C  iai 


Introduction 

In  the  grading  of  timber  for  beams 
it  is  useful  to  have  a  quantitative 
measure  of  the  effect  of  end  splits  on. 
the  strength  of  such  beams.  The 
method  presented  gives  a  way  of 
predicting  when  an  end  crack  will 
propagate  resulting  in  a  loss  of 
stiffness  and  strength.  Linear  elastic 
fracture  mechanics  has  been  applied 
successfully  to  predict  brittle  failure 
(crack  propagation)  in  wood  (2.8.  li)’ 
Although  the  title  problem  has  been 
examined  (2),  this  report  uses  a 
different  method  of  solution  and 
simplifies  the  results  while  investigat¬ 
ing  the  effect  of  different  orthotropic 
parameters. 

Before  proceeding,  linear  elastic 
fracture  mechanics  will  be  over¬ 
viewed.  The  propagation  of  a  planar 
crack  can  be  described,  in  the  most 
general  case,  by  the  combination  of 


three  independent  modes  as  shown 
in  figure  1.  The  first  mode,  or 
opening  mode,  has  the  crack  sur¬ 
faces  move  in  a  normal  direction  with 
respect  to  the  surface  planes.  The 
second,  sliding,  or  forward  shear, 
mode  has  the  surfaces  move  in  their 
respective  planes  in  the  direction 
parallel  to  collinear  crack  propaga¬ 
tion.  The  third,  tearing,  or  transverse 
shear,  mode  has  the  surfaces  moving 
in  their  respective  planes  but 
perpendicular  to  collinear  crack 
propagation.  The  opening  mode 
results  from  inplane  (a  plane 
perpendicular  to  the  crack  front) 
symmetric  (symmetric  about  the 
plane  between  crack  plane  surfaces) 
loading  while  the  sliding  mode 
results  from  inplane  skew-symmetric 
loading.  The  third  mode  results  from 
out  of  or  antiplane  loading  conditions 
(loading  parallel  to  crack  front). 

The  theoretical  stresses  (see  fig.  2 


for  positive  stresses)  for  a  small 
region  directly  ahead  of  the  crack 
(the  crack  surfaces  coincident  with 
the  negative  x,  z  half  plane)  for  the 
three  modes  are  reduced  from  Sih  et 
al.  (15). 


Mode  I 


Mode  II 


<>,  =  K|i)i»i2(2jrx) <i,  =  0 

nv  =  K|(2 7rx)  "*  nv  =  0 

■«v  =  0  rxv  =  Ku(2ttX) 

Mode  III 

r„,  =  0 

V  =  Kin(2irx) 

(>;  =  0 


'  Maintained  at  Madison.  Wts  .  in  cooperation 
with  the  University  ot  Wisconsin. 

’  Underlined  numbers  In  parentheses  refer  to 
literature  cited  at  ehd  ot  report 


Figure  1.— Three  basic  modes  of  crack  surface  displacement. 
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MODE  m 
"TEARING" 


<§>  a  -  Va  #  R 
<g>  a  -■  Maz  Ra 


Figure  2.— Two-dimensional  coordi¬ 
nate  system  around  a  crack  tip 
with  positive  stresses. 
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Figure  3.— End-split  beam  under  concentrated  load  (superposition  of  two 
basic  cases). 
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where  are  orthotropic  param¬ 
eters,  that  is,  the  13' s  characterize  the 
degree  of  orthotropy  and  are  equal  to 
one  for  isotropy  (see  appendix  A), 
and  the  K’s  are  stress  intensity  fac¬ 
tors  (SIF’s)  in  I,  II,  or  III  loading 
modes.1  The  SIF’s  are  linearly 
dependent  on  remote  boundary 
conditions,  strongly  dependent  on 
geometry,  and  weakly  dependent  on 
orthotropic  parameters  for  finite 
bodies  (the  dependence  on  ortho¬ 
tropic  parameters  becomes  stronger 
as  the  distance  ahead  of  the  crack 
and  a  free  surface  becomes  smaller). 
In  single-mode  loading  and  collinear 
crack  propagation,  the  crack 
propagates  when  the  SIF  reaches  a 
critical  value  (K,c,  K„c,  or  K,„c). 

Theoretical  Considerations 

A  boundary  value  collocation 
method  was  employed  to  solve  for 
the  theoretical  SIF  K„.  First  a  series 
solution  was  formulated  which  satis¬ 
fies  identically  the  field  equation  (a 
generalization  of  the  biharmonic 
equation)  and  the  stress-free  condi¬ 
tion  on  the  crack  surfaces  (see 
appendix  A).  The  unknown  coeffi¬ 
cients  of  the  series  solution  were 
solved  using  boundary  value  colloca¬ 
tion  (see  appendix  B).  Boundary 
values  were  resultant  forces  and 
moment  (see  eq.  A5,  A6).  This 
method  was  employed  because  (1) 
after  the  coefficients  are  obtained, 


stresses  and  strains  can  be  calcu¬ 
lated  anywhere  In  the  body,  (2)  the 
first  coefficient  is  linearly  related  to 
the  SIF,  (3)  the  set  of  linear  equations 
result  in  a  dense  matrix  requiring 
minimal  computer  storage  space,  and 
(4)  more  than  one  problem  (different 
boundary  values)  can  be  solved 
simultaneously.  The  number  of 
boundary  points  was  increased  until 
the  SIF  (first  coefficient)  converged, 
(e.g.,  in  the  concentrated  load  case, 
increasing  the  number  of  boundary 
points  (each  with  three  stress  condi¬ 
tions)  from  24  to  32  changed  the  SIF 


less  than  0.8  pet).  A  more  detailed 
description  of  the  theoretical  analysis 
is  given  in  appendices  A  and  B. 

Theoretical  Problem 
and  Results 

End-Split  Beam  Under 

A  Concentrated  Load 

The  problem  considered  is  an  end- 
split  beam  (end  split  horizontal,  full 
width,  and  at  middepth)  as  shown  in 
figure  3.  Assume  no  tractions  on  the 

'  Some  authors  us#  ki.  kj,  snd  kj  as  SIF's 
wtisr*  K|  *  k)  t  ,  K||  a  kj  v  *  ,  snd 
k«i  =  kj  . 
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only  one  more  subcase  is  added  than 
the  concentrated  load  case.  The  K„ 
for  a  uniform  load  can  then  be  esti¬ 
mated  (conservatively)  as: 


K„ 


A  — M_  +  B  V— 
bd’'’.  bdwl 


(2) 


+  C_* 
bd  ’ 


where  w  is  the  uniform  loading  and  C 
is  a  third  constant  dependent  on  the 
orthotropic  parameters.  Note  equa¬ 
tion  (2)  reduces  to  equation  (1)  for  the 
concentrated  load  case  (e  g.,  w  =  0). 
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Effect  of  Different  Orthotropic 
Parameters 
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Figure  4. — End-split  beam  under  uniform  load  (superposition  of  four  basic 
cases). 
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surfaces  of  the  split.  The  problem 
can  be  solved  by  superposing  the 
solutions  of  two  problems  shown  in 
figure  3.  The  two  cases  are  solved 
and  then  the  SIF  can  be  written 
(using  superposition): 


K„ 


A-M-  +  B_V_ 
bd*"  bd1" 


(D 


where  the  internal  moment  and  shear, 
M  and  V,  are  calculated  at  the  crack 
tip  using  elementary  methods,  A  and 
B  are  constants  dependent  on  the 


orthotropic  parameters,  and  b  and  d 
are  the  beam  width  and  depth, 
respectively. 

End-Split  Beam  Under 

A  Uniform  Load 

This  problem  is  shown  in  figure  4 
with  subsequent  breakdowns  (again 
assume  no  crack  surface  tractions). 
The  subcase  with  symmetric  loading 
gives  a  negative  opening  mode  SIF 
which  would  produce  frictional  forces 
between  the  two  crack  planes:  there¬ 
fore,  as  a  conservative  estimate,  K, 
will  be  taken  as  zero.  As  can  be  seen. 


Orthotropic  parameters  for 
different  wood  species  and  grain 
orientations  and  plane  stress  condi¬ 
tions  are  shown  in  figure  5  (see  eq. 
(A2),  (A3)).  The  coefficients  A.  B,  and 
C  (see  eq.  (2))  have  been  calculated 
for  four  extreme  and  one  average 
combination  of  /S,  and  /J2  shown  in 
figure  5.  The  coefficients  are  given  in 
table  1  and  their  respective  concen¬ 
trated  load  equations  (eq.  (1))  are 
plotted  in  figure  6.  An  equation  to  be 
used  for  any  wood  species  and  any 
grain  orientation  (constants  have  two 
significant  figures)  under  concen¬ 
trated  or  uniform  load  conditions  is 


K„  = 


-2.8 


M 

bd’" 


0.72 


V 

bd”' 


+  0.06  -Jr.  ...  (3) 

bd  ’'J 


Equation  (3)  is  good  (within  11  pet) 
for  all  combinations  in  figure  6. 
Combined  loading  conditions  can  be 
considered  by  summing  up  the  Kn's 
for  each  loading  condition  yielding  a 
total  Km.  Failure  occurs  when  ab¬ 
solute  total  K„  reaches  a  critical 
value,  Kllc.  Therefore  equation  (3)  with 
an  appropriate  K„c  becomes  a  single 
failure  criteria  for  end-split  beams. 

End-Split  Beams  Concentrated 

Load  Experiments 

The  experimental  data  are  from  an 
unpublished  report  by  C.  B.  Norris 
and  E.  C.  O.  Erickson,  done  in  1951. 
The  following  section  is  verbatim* 
from  that  report: 


Description  of  Tests 

“Two  groups  of  substantially  clear 
Douglas-fir  beams  (7.5-  by  1.62-in 

‘  Figures,  tables,  and  references  are  renum 
bered 


3 


ft 


Table  1.— Rang*  of  eoaffleianta  In  an  aquation  of  KM  as  a  function  of  loading 
conditions. 
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ii 

~w 

+  B  v 

+  r.  w 

bd*'1 

bd’" 

bd-"1 

Orthotropic  parameters 

Coefficients 

0, 

@2 

A 

B 

C 

0.75 

3.0 

-3.112 

-0.657 

0.057 

1.5 

3.0 

-2.864 

-  .674 

.060 

1.5 

4.25 

-  2.577 

-  .779 

.060 

.75 

4.25 

-  2.728 

-  .769 

.068 

1.125 

3.625 

-2.790 

-  .723 

.062 

25 

20 

15 
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Figure  6. — Graph  of  dimensionless  KM  as  a  function  of  relative  crack  length 
for  an  end-split  beam  under  concentrated  load. 
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Figure  5.— Orthotropic  parameters  for 
different  wood  species  for  the  TL, 
RL  systems  of  crack  propagation. 
(First  letter  refers  to  normal  to 
crack  planes  or  the  y  direction;  the 
second  letter  refers  to  direction  of 
crack  propagation  or  the  x  direc¬ 
tion.)  Computed  from  elastic  con¬ 
stants  in  (j)  or  ‘obtained  from  J.  0. 
Barrett  private  communication. 
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cross  section)  were  tested  as  shown 
in  figure  7.  Saw  kerfs  about  1/8-inch 
wide  and  of  various  lengths  (see  table 
2)  were  cut  at  the  centerline  of  the 
beams  as  shown  in  the  figure.  These 
kerfs  were  filled  with  smooth  splines 
of  wood  which  were  well  paraffined 
to  reduce  friction.  The  beams  were 
supported  on  rockers  and  protected 
from  indentation  of  the  knife  edges 
by  steel  plates  1-inch  thick  and 
extending  3  inches  each  side  of  the 
reaction  in  the  direction  of  the  length 
of  the  beams.  The  beams  were 
laterally  supported,  near  their  ends, 
during  test  to  keep  them  in  a  vertical 
plane.  These  supports  did  not 
hamper  the  deflection  of  the  beams. 

“Pillow  blocks  of  hardwood  were 
used  to  transfer  the  loads  to  the 
beams.  For  the  beams  of  group  1,  the 
pillow  blocks  were  out  to  a  24-inch 
radius;  and  for  the  beams  of  group  2, 
to  an  18-inch  radius.  The  speed  of  the 
testing  machine  was  0.194  inch  per 
minute  for  group  1  and  0.096  inch  per 
minute  for  group  2.  The  deflections 
were  measured  to  the  nearest  0.01 
inch  at  load  increments  varying  from 
200  to  500  pounds. 


"After  the  beams  were  tested,  two 
block-shear  specimens  were  cut  from 
each  of  them  in  the  neighborhood  of 
the  end  of  the  saw  kerfs  and,  subse¬ 
quently,  tested  Q). 

Presentation  of  Data 

"The  data  obtained  from  the  tests 
and  some  calculated  values  are  given 
in  table  2.  The  first  column  contains 
the  beam  number.  Beams  3,  9,  and  1C 
were  omitted  from  the  table  because 
they  were  cross  grained.  Beam  15 
was  omitted  because  the  shear 
failure  was  so  gradual  that  it  was  not 
noted  until  the  specimen  was 
removed  from  the  testing  machine. 
The  second  column  gives  the  length 
of  the  saw  kerf  measured  from  the 


position  of  the  reaction  toward  the 
center  of  the  beam,  as  shown  in 
figure  7.  The  third  column  gives  the 
distance  from  the  end  of  the  kerf  to 
the  position  of  the  load.  Column  4 
gives  the  largest  loads  the  beams 
withstood;  and  column  5  the  first 
maximum  load.  These  two  loads  are 
identical  except  for  four  beams.  In 
these  four  tests  the  shear  failures 
were  not  completely  destructive;  the 
beams  failed  and  the  load  dropped, 
and,  as  the  test  proceeded,  the  load 
increased  to  a  value  greater  than  the 
first  maximum. 

"Column  6  of  the  table  lists  the 
value  of  the  reaction  at  the  first 
maximum  load,  and  column  7  the 
shear  stress  computed  by  the 
formula 
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Figure  7  — Sketches  of  beams  tested. 
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into  account  the  depth  of  the  beam 
It  would  be  expected  that  the 
presence  of  this  kerf  would  not 
weaken  the  beam 

"The  saw  kerfs  in  beams  17  and  18 
extended  to  within  7  and  3  inches  ol 
the  position  of  the  load  These  two 
beams  are  not  expected  to  satisfy  the 
suggested  equation  because  the 
ends  of  the  kerfs  are  too  close  to  the 
position  of  the  load." 

Comparison  of  Experiment 

with  Theory 

To  analyze  the  experimental  data,  a 
theoretical  equation  similar  to  equa¬ 
tion  (3)  is  needed.  (Eq.  (3)  is  a 
suggested  design  equation,  too 
general  when  comparing  theory  with 
experiment.)  Using  the  orthotropic 
constants  published  (see  (5))  for 
Douglas-  fir  (radial-longitudinal,  RL) 
and  a  stress  function  boundary  value 
collocation  method  (see  appendices 
A  and  8).  the  equation  for  the  mode  II 
SIF  was  found  to  be 


Table  2.— Pertinent  data  from  the  tests  and  some  calculated  values 


Speci¬ 

men 

No. 

a 

c-a 

Maxi¬ 

mum 

load 

First 

maxi¬ 

mum 

load 

Reac¬ 

tion 

rc 

Results  of  block-shear 
tests 

Speci¬ 
men 
No.  1 

Speci¬ 

men 

No.  2 

Aver¬ 

age 

(D 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

In. 

Lb. 

Psi 

GROUP  1 

5 

0 

36 

8.420 

8,420 

4,210 

516 

1,480 

1,472 

1,476 

1 

4 

32 

8,150 

8,150 

4,075 

500 

1,416 

1,421 

1,418 

6 

8 

28 

4,700 

4,700 

2,350 

288 

1,327 

1,225 

1,276 

2 

12 

24 

3,690 

3,340 

1,670 

205 

1,342 

1,319 

1,330 

4 

12 

24 

2,670 

2,670 

1,335 

164 

1,171 

1.140 

1,156 

7 

16 

20 

2,630 

2,630 

1,315 

161 

1,389 

1,616 

1,502 

8 

20 

16 

2,510 

1,550 

775 

95 

1,336 

1,181 

1,258 

GROUP 2 

11 

5 

19 

7,170 

7,170 

4,710 

578 

1,519 

1,526 

1,522 

13 

5 

19 

5,240 

5,240 

3,440 

428 

1,435 

1,302 

1,368 

16 

9 

15 

3,225 

3,225 

2,120 

260 

1,352 

1,222 

1,287 

19 

9 

15 

3,090 

3,090 

2,030 

249 

1,277 

1,130 

1,204 

12 

13 

11 

3,000 

2,500 

1,640 

201 

1,339 

1,363 

1,351 

’14 

13 

11 

2,360 

2,360 

1,550 

195 

1,233 

1,187 

1,210 

17 

17 

7 

3,900 

2,700 

1,770 

217 

1,326 

1,470 

1,398 

18 

21 

3 

3,600 

3,600 

2,360 

290 

1,181 

1,216 

1.198 

•  Specimen  14  had  a  depth  ol  7.33  In. 


Columns  8,  9,  and  10  give  the  results 
of  the  block-shear  tests  made  on 
specimens  cut  from  the  beam  in  the 
general  location  of  the  end  of  the 
saw  kerf. 

"The  saw  kerf  in  beam  5  was  cut 


through  the  4-inch  overhang  of  the 
beam  just  to  the  position  of  the 
reaction.  This  kerf  did  not  cause  a 
shear  failure.  The  beam  failed  at  a 
modulus  of  rupture  of  9,930  pounds 
per  square  inch,  which  is  consistent 
with  the  average  value  of  11.200  given 
for  coast-type  Douglas-flr  (10),  taking 


K„  =  -2.785  M  _  0  731  V 
bd’  >  bd’  ' 

+  0.0646  W  (4) 

bd  ■'» 

For  a  concentrated  load  case, 
equation  (4)  can  be  rewritten  to 

Kh  =  [-2.785(4)  -  0.731]  R  (5) 
d  bd1  ■ 

where 

K|,  is  the  sliding  mode  SIF 
a  is  crack  length  (column  2, 
table  2) 

d  is  beam  depth  (7.5  in.) 
b  is  beam  width  (1.62  in.) 

R  is  reaction  at  support 
(column  6.  table  2). 

Using  the  experimental  data  and  the 
theoretical  equation  (5).  values  of  KMc 
were  calculated  for  12  beams  (beams 
5.  17.  and  18  were  not  considered 
because  of  reasons  noted  earlier)  and 
are  listed  in  table  Cl.  Because  group 
1  beams  had  two  end  splits  and 
group  2  beams  had  only  one,  some 
adjustment  between  the  two  groups 
had  to  be  made  to  obtain  an  average 
critical  SIF  for  a  beam  with  only 
one  end  split.  This  adjustment  is 
explained  in  appendix  C.  From  the 
data  of  the  12  beams  and  statistical 
calculations,  KJ,C  was  found  to  be 
2,035  psi  x  in.  with  a 
coefficient  of  variation  (COV)  =  15 
percent,  for  a  beam  with  one  end 
sp'it.  This  compares  to  the  value 
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obtained  by  Barrett  and  Foschi  (2)  tor 
western  hemlock  of  2,107  psi 

v  in. .  Figure  8  graphs  -1  _B_ 

2  bd 


versus  A  using  equation  (5) 
d 


(w  =  0}  and  =  2,035  and  1,885 
psi  v  in.  for  single  and  double 
end-split  beams,  respectively.  Also 
plotted  in  figure  8  are  the  experi¬ 
mental  data  showing  the  excellent 
agreement  between  experiment  and 
theory.  If  the  generalized  species 
equation  (3)  is  used  instead  of 
equation  (4)  with  the  test  data,  the 
would  be  calculated  to  be  about 
Vi  percent  higher. 


Conclusions 


1.  Boundary  value  collocation  can 
be  used  to  accurately  calculate  the 
shear  mode  SIF  in  orthotropic  wood 
beams. 


2.  Using  superposition  principles, 
one  single  equation  can  be  written  to 
calculate  K„  for  an  end-split  beam  for 
either  concentrated  or  uniform  load 
cases. 

3.  Complex  loadings  can  be 
handled  by  superposition  of  Ku’s. 

4.  The  effect  of  ditterent 
orthotropic  parameters  on  the 
simplified  equation  is  small  enough 
so  that  one  equation  can  be  used  for 
a  majority  of  structural  wood  species. 

5.  Using  the  BVC  method  and 
species  specific  coefficients  in  the 
simplified  equation.  for  a  species 
can  be  determined  experimentally. 

6.  Douglas-fir  KL  In  the  RL  system 
was  determined  to  be  2,035  psi  \  in. 
at  a  moisture  content  of  11  percent. 

7.  A  proposed  simple  failure  equa¬ 
tion  for  end-split  beams  is 


Kiic  s 
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Figure  8.— Comparison  of  experimental  data  with  theory,  for  end-split  beams 
under  concentrated  loads  (•—single-end  split,  o-  -double-end  split). 
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splits.  A  simple  failure  equation  for  end-split  beams  splits,  A  simple  failure  equation  for  end-split  beams 

Is  also  proposed.  is  also  proposed. 
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Appendix  A 

Formulation  o<  tho  Compton  Sort**  Stroos  Functions  tor  Crocked  Orthotropic  Bodies 

Sih  et  al.  (IS)  have  given  a  general  complex  series  stress  (unction  (or  a  cracked  rectiiinearly  anisotropic 
body  These  (unctions  simplify  when  the  crack  is  parallel  to  a  material  axis  Leading  up  to  these  functions 
consider  first  an  orthotrop  c  body  with  a  crack  coincident  with  the  negative  x-axis  (tig  Al)  and  a  material  axis 
The  equilibrium  equations  (no  body  forces): 


-*’>  =  0 

z: 

<1x 

<fy 

,iy  ,lx 

and  compatibility  equation: 

•V'v  +  >’S  = 

■fy!  <ty' 

•V,M 

3x«iy 

are  combined  into  one  field  equation: 

a« 

.VU 

,tx* 

+  (2a, :  +  a**)  (1^yJ 

*  %  -  0 

using  the  following  relations 

«\ 

,VU  ,  .VU 

3y'  l’‘  ,ix»  T" 

,i;U 

,*x.»y 

—  am',  +  a,.,\ 


and 


*  xv  “ 

The  coefficients  aM  are  related  to  the  engineering  constants  by 


Ev  •  =  C  '  "  l  '  **  =  g[ 


an  -  £  .  at:  — 


for  plane  stress  and 


a,,  =  1  ,  a,..  =  '*»'  " 

tx 


=  1  Ya'‘«  a«  = 


1 


(All 


(A2) 


for  plane  strain 

Specializing  the  rectiiinearly  anisotropic  case  of  (3.9.1J.13).  the  solution  to  the  orthotropic  field  equation  (An 
in  terms  of  analytic  functions  F,of  complex  variables  z,  is 

U(z,,Zi)  =  2  Re[F,(z,l  +  F..(z..)) 

where 


z,  =  x  +  id,y  =  r(cos  t>  +  id,  sin  t>\  -  x  -s  r*  ^  x 


where 


i  =  x  -  1 


Re  means  the  real  part  of,  and  the  d,  s  satisfy  the  characteristic  equation 


Introducing  the  notation 


~  (2St?  +  +  a;j  -  0. 


(A3) 
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The  remaining  coefficients  (a,„'s  and  btn's)  of  the  complex  series  stress  function  are  obtained  using  a  boundary 
value  collocation  method  as  described  in  appendix  B. 

The  leading  coefficient  X11  is  related  to  the  SIF's  Kh  Kt|  by  the  relationship: 

K,  =  -  a„ 

and 

K|,  =  b„  -JBir  (0,  -  02). 


Appendix  B 

Boundary  Value  Collocation  and  Numerical  Method 


In  a  boundary  value  collocation  method  (or  point  matching)  the  form  of  the  solution  satisfies  the  governing 
field  equation  everywhere  in  the  interior  of  the  body.  The  constants  or  coefficients  appearing  in  the  solution  are 
obtained  by  evaluating  the  solution  at  a  finite  number  of  boundary  points.  Infinite  series  solutions  are  first 
truncated  such  that  the  number  of  unknown  constants  are  at  most  equal  to  the  number  of  boundary  conditions 
to  be  evaluated.  If  more  boundary  conditions  are  specified  than  unknown  coefficients,  the  unknowns  can  be 
solved  for  in  a  least-squares  sense.  Solving  in  a  least-squares  sense  has  the  advantage  of  reducing  off  point 
errors,  though  making  the  linear  system  ot  equations  more  ill-conditioned. 

To  apply  the  boundary  value  collocation  method  to  a  cracked  orthotropic  body,  the  complex  stress  function 
(see  appendix  A)  and  associated  integrals  are  approximated  by: 

M  n  -  i 

M  Xi 

4>,'<z,)  s  E  ,  A"  W 

n-i  j 


and 


M 

F'<zi>  *?„r 


(Zi)* 


\n _ 

'<"  *  1) 
i  2 

These  functions  are  now  used  (according  to  eq.  (A4)  to  (A6)  to  specify  stress,  displacement,  or  resultants  at  a 
discrete  number  of  boundary  points,  and  the  resulting  linear  equations  can  be  written  in  matrix  notation  as 

A  X  =  B 


where  X  contains  M  (not  to  be  confused  with  moment  M  in  text  and  appendix  A)  unknown  complex  constants  \|n. 
matrix  A  contains  functions  of  position  and  0l’s,  and  B  contains  N  corresponding  boundary  stress,  displacement, 
or  resultant  values.  For  the  orthotropic  crack-alined  case  the  Re[\nJ  =  a,„  corresponds  to  symmetric  (about  the  x 
axis)  in-plane  loading  only,  while  lm[X|n]  =  b,„  corresponds  to  skew-symmetric  in-plane  loading.  Thus  the  general 
case  is  really  the  sum  of  two  separate  problems: 

A  Re(\]  =  Bsym  and  A  lm[\]  =  BSKEW 

which  are  solved  separately.  Each  case  automatically  satisfies  appropriate  stress  derivatives  ahead  of  the  crack 
and  because  of  geometric  symmetry  only  the  region  y  >  0  has  to  be  considered.  The  procedure  used  will  be 
given  and  then  the  rationale  for  each  procedural  step  will  be  explained  later. 

Begin  first  with  an  overspecified  system 

AX  =  B, 

where 

A  is  an  N  X  M  known  matrix 
X  is  an  M  X  1  unknown  column  vector 

and 

B  is  an  N  X  1  known  column  vector 

and 


M  s  N. 

First  postmultlpiy  A  by  an  (M  x  M)  diagonal  matrix  0  to  get 

AD  D  '  X  =  ADy  =  Cy  =  B  (B1) 
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and  premultiply  C  and  B  by  the  transpose  of  C  knowing  that 

CTCy  =  CTB. 

The  square  (M  x  M)  matrix  CrC  is  equal  to  TTT  where  T  is  an  upper  triangular  (M  x  M)  and  is  obtained  by 
triangular  reduction  by  an  orthonormal  matrix  R  such  that 

RA  =  (l) 

(see  (12)).  The  matrix  equation  is  now 

rTy  =  CTB 

Taking  the  inverse  of  the  left  side  gives  the  matrix  equation  for  the  unknown  column  vector 

X  =  DT’TTCTB 

Both  the  matrices  X  and  B  can  be  generalized  to  K  columns  so  that  a  K  number  of  problems  (with  the  same 
orthotropy  and  geometry  and  different  loadings)  can  be  solved  simultaneously. 

The  matrix  A  is  scaled  by  postmultiplying  it  by  the  diagonal  matrix  0  in  equation  (B1)  to  achieve  greater 
accuracy  in  the  solution  than  when  it  is  not  scaled.  The  resulting  scaled  matrix  is  required  to  have  column 
Euclidean  norms  approximately  equal  to  unity  and  the  diagonal  elements  of  0  are  powers  of  two  to  avoid 
introducing  computer  roundoff  errors.  Finding  T  from  A  directly  by  triangular  reduction  by  an  orthonormal  matrix 
is  believed  to  be  numerically  more  stable  than  by  Gram-Schmidt  triangular  reduction  (12)  and  requires  less 
computation  than  using  Cholesky  factorization  of  CTC  into  TtT.  Finally  then,  the  only  matrix  inversion  is  the 
simple  inversion  of  a  triangular  matrix  (T).  It  was  found  that  scaling  and  triangular  reduction  increased  the 
accuracy  of  X  and  in  some  instances  (M  large)  gave  answers  when  elimination  methods,  double  precision,  and 
iteration  improvement  did  not. 

Appendix  C 

Statistics 

Using  the  experimental  data  from  table  2  and  the  theoretical  equation  (4).  table  Cl  is  constructed 
Since  in  group  1  the  weaker  of  two  K„  was  measured,  the  data  cannot  be  lumped  together  without  some 
algebraic  adjustment  first.  Assume  that  for  a  single  set  of  conditions  (one  split,  etc.)  K„c  is  normally  distributed 
N<mi.  <>,'):  assume  also  that  for  two  splits  the  measured  Knc  is  the  lower  value  of  two  independent  K„c,  each 
having  the  same  distribution;  then,  the  probability  density  function  (p.d.f.)  of  this  measured  KMl.  is  (using  series 
reliability  (4)) 

9?<r.  Mi.  «'i>  =  2[1-FN(y;  mi.  <>,))  fN(y:  mi.  "i)  (Cl] 


mi.  <’i)  — 


Ot\  2  IT 


Fn(x;  /I,, ..,)  =  |^fN(y;  o,)dy 

Using  the  method  of  maximum  likelihood  (6),  we  write  a  likelihood  function  L  as 

L0*i.  0i.'  x,.  Xj . x8,  y1t  y2 . y6> 

0 

=  j  i  9 ?(y •  i  in.  °i)  *n(x,;  mi.  0|) 

where  x,  are  observed  values  from  group  2  (column  3.  table  Cl)  and  y,  are  observed  values  from  group  1  (column 
4,  table  Cl).  The  maximum  likelihood  estimates  for  and  «,  are  found  such  that 


!?l  =  0  and  ‘K  =  0, 

Ofl  f  (Ml  J  i 

This  was  done  with  the  data  from  table  Cl  resulting  In 

Mi  =  2,035  psi  v  in. 

and 

fli  =  301  psl  vMn.  • 

These  are  estimates  of  the  mean  and  standard  deviation  for  a  single-split  beam. 
The  expected  value  for  the  group  1  set  is 

m?  =  i  -  x  g*(x)dx 


With  a  change  of  variables  equation  (C2)  becomes 

»■}  =  r.(x.<,  +  («,)  2(1  -  <t>(x)]  o<x)dx 

where  o£x)  is  the  p.d.f.  and  *(x)  is  the  c.d.f.  (cumulative  distribution  function)  of  the  standardized  normal 
distribution.  Equation  (C3)  can  be  evaluated: 

><2  =  ■■,!*  x2[1  - <t>(x)|  o(x)dx  +  ,.,|1  -♦(x))’ 

=  Co,  +  /‘i 


The  variance  is 


but 


of  =  (1  x'gj(x)dx  -  ii 


("  x,g,(x)dx  =  (1  (x«>,  +  /<,)'  2(1  -  <t*(x)|  o<x)dx 
=  of  It*  x'<Mx)dx  +  x*2(Vi  -  <h(x))tf>(x)dx| 

+  2(1,0, (*  x2(1  -  «t>(x)]  o(x)dx  +  2|1  -  <t>(x)]o(x)dx 

=  <»;|i  +  o]  +  2(i, o,c  +  (i? 

resulting  in 

of  =  of  +  2(1, o,C  +  (if  —  (if 


Using  equation  (C4)  gives 


of  =  |1  -C1)..; 


The  constant  C  was  calculated  by  numerical  integration  to  be  -0.5642  yielding: 


ft,  =  (i,  -  0.5642iT  =  1,865  psi  \  in. 


and 


(C3) 


(C4) 


(C5) 


o,  =  0.8256o,  =  249  psi  \  in.  (C6) 

Therefore,  if  a  beam  has  one  split  which  would  always  fail  first,  (i,  and  o,  would  describe  the  mean  and 
variance  of  K(l<„  and  if  a  beam  has  two  splits  (each  with  the  same  Kn),  ,i:  and  o.,  would  describe  the  mean  and 
variance  of  the  observed  K„r. 


Table  Cl.— Computed  Kllc  trom  experiments 


Specimen 

number 

K||(„  x. 

Kn,..  y, 

Specimen 

number 

Kn,-  X,  Kn,..  y, 

1 

GROUP  1 

2.036 

11 

GROUP 2 

2.747  - 

6 

— 

1.961 

13 

2,006  - 

2 

— 

1.952 

16 

1.946  - 

4 

— 

1,561 

19 

1.864  - 

7 

— 

1.978 

12 

2.055  - 

8 

— 

1,425 

14 

2,004  - 

2.0-12/79 
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